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Abstract 

We derive dual representations for 0(N) and CP(N-l) models on the lattice. In terms of the dual 
variables the partition sums have only real and positive contributions also at finite chemical potential. 
Thus the complex action problem of the conventional formulation is overcome and using the dual 
variables Monte Carlo simulations are possible at arbitrary chemical potential. 
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1 Introductory remarks 

A very visible shortcoming of lattice QCD is the current inability to properly deal with QCD at finite 
density. This is not a fundamental conceptual problem, but more of a technical challenge: At finite 
chemical potential fj, the action S is complex and the Euclidean Boltzmann factor e~^ does not have 
a probability interpretation which is necessary for a Monte Carlo simulation. As a matter of fact, 
the complex action problem appears in most field theories at finite chemical potential, both on the 
lattice and in the continuum, and for fermionic as well as for bosonic degrees of freedom. In these 
systems importance sampling methods are thus not able to identify ‘relevant configurations' among the 
fundamental degrees of freedom, which might hint at more effective degrees of freedom allowing for a 
better understanding of physical effects at finite density. As a side remark we also note that the addition 
of a topological term leads to a complex action problem, similar to the one at finite //. 

A very powerful approach to completely solving the complex action problem is to exactly rewrite the 
partition sum in terms of new degrees of freedom such that the partition sum has only real and positive 
contributions. These new degrees of freedom are often referred to as dual variables which are given by 
loops for matter fields and by surfaces for gauge fields. In recent years several interesting lattice field 
theories at finite chemical potential or with a topological term were rewritten in terms of dual variables 
(for two reviews see, e.g., [Il[2])- In the dual formulation the system can, e.g., be updated efficiently 
with worm algorithms of the Prokofev-Svistunov type 0, which also has a generalization to the surfaces 
that dualize the gauge fields [4]. 
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An important lesson that has been learned when developing the dual approach is that there is no 
general strategy for finding a dual representation. Every type of model has to be analyzed individually 
and in particular the type of symmetry plays an important role since it directly determines the constraints 
for the dual variables which in turn determine the structure of the dual loops and surfaces. A second 
insight that has emerged is that a dual representation is not unique - there might be several different dual 
representations for the same lattice field theory. The different dual representations may have different 
properties, in particular with respect to a possible dual Monte Carlo simulation. 

An important class of systems are lattice field theories with O(A^) symmetric actions and the closely 
related CP(A^ —1) models. These models are especially interesting in two dimensions because of their 
close resemblance to gauge theories. In particular both 0{N) and CP(A^ —1) in two (1+1) dimensions 
are asymptotically free and have a dynamically generated mass gap proportional to the strong scale of 
the theory. In addition 0(3) and CP(A^—1) for any N have topological charge and instantons and allow 
for a 6 term, making them even more attractive for understanding the gauge dynamics. 

We mention in passing that various chemical potentials which we discuss here have a close similarity 
with twisted boundary conditions for the fields of the model (the chemical potentials can be analytically 
continued to the twists), and were studied extensively in the context of fractional topological charge and 
resurgence [5HI5]. 

For models with these symmetries, dual variants were discussed in the literature [Ti riBl - [22] . Here we 
revisit the problem of dualizing 0{N) and CP(A^ —1) models with finite chemical potential. We show 
that both systems allow for a dualization along the same lines using N sets of dual flux variables for the 
0{N) case and 2N sets for CP(A^—1). It is interesting to note that the dualizations presented here differ 
fundamentally from other suggestions (which only considered the case without chemical potential): For 
the O(A^) model a dualization with only one set of dual variables was presented in [17] (we use N sets of 
dual variables), while for the CP(A^—1) case variants with 2 [18] and with [I] sets of dual variables 
were discussed (we use 2N sets). 

The fact that rather different exact dual representations are available for the same model is interesting 
for at least two reasons: Identifying a real and positive dual representation is only the first step of the dual 
approach. One also has to find a suitable Monte Carlo algorithm for updating the dual degrees of freedom 
which are subject to constraints that have to be obeyed during the update. Different representations 
that use completely different dual variables will necessarily give rise to Monte Carlo update schemes 
with different properties and different efficiency. Exploring the space of possible dual representations is 
certainly important for identifying the most efficient schemes. The second reason for exploring different 
dual representations is that currently for the probably most interesting case of non-abelian gauge systems 
so far only dual representations that also have negative weights were found [23H3Q] . However, with every 
new dualization strategy we enlarge our toolbox and get closer to the goal of finding real and positive 
representations suitable for dual Monte Carlo simulations. 

2 Lattice formulation of 0{N) and CP{N—1) models with a chem¬ 
ical potential 

In this section we first discuss the conventional lattice discretization for the O(A^) model with chemical 
potential and then for the CP(A^—1) model. For vanishing chemical potential the action of the nonlinear 
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0{N) model on the lattice is given by 


d 

5 ’[^] = -J'^'^r{x)-r{x + v) , 

xgA v=l 


( 1 ) 


where the first sum runs over all sites of the d-dimensional lattice A with size V = x Nt and 

periodic boundary conditions in all directions. The second sum runs from u = 1 to v = d, and by 
i> we denote the unit vector in r/-direction. The direction v = d will be referred to as the Euclidean 
time direction. The dynamical degrees of freedom are real A^-component spin vectors r{x) assigned to 
the sites of the lattice. The vectors obey the constraints f(x)^ = 1, i.e., they are constrained to the 
N — I sphere S^~^. The action ([1]) consists of a ferromagnetic nearest neighbor interaction of the spin 
vectors with a coupling strength J (the lattice spacing a is set to a = 1 throughout this paper). The 
partition sum is given by Z = / P[r ] exp(—S'[r ]) and the coupling J is positive (ferromagnetic case). 
The integration measure V[r\ is the product of the 0{N) invariant measures over at each site 

(see ([6]) below). 

Due to the global O(A^) symmetry the model has conserved charges which can be coupled to chemical 
potentials. Following [2Tl|3l] we introduce a chemical potential /i to one of them (all other cases can 
be obtained by adding several such "fundamental” chemical potentials [2I])- The chemical potential 
gives a different weight to forward and backward nearest neighbor terms in time direction {v = d). The 
action with chemical potential thus reads 

d 

5[f] = ^ [rj_(x) • rj_(x + z)) 

xsA i/=l 


+ 


■ ri 2 {x) ruix + v) 


( 2 ) 


We use a notation where the first two components are combined into a complex number ri{x) + ir 2 {x) = 
ri 2 (x) and r_\_{x) contains the N — 2 components which do not couple to the chemical potential, 


r± = (rsix), r4(x),.. .rjv(x)) . (3) 

We require r’i 2 (a:)^ + f±(x)^ = 1 and thus for /i = 0 the action ([2]) reduces to ([T]). We stress that in 
the conventional form the action is complex for /x / 0 and a Monte Carlo simulation suffers from the 
complex action problem. 


Let us now come to the CP(A^—1) model where various different lattice discretizations were discussed 
(see, e.g., [I8ll32] ). We here use a variant where the interaction is introduced with an auxiliary link field 
Ui,{x). For CP(A^ —1) the chemical potentials can be added for the charges associated with the global 
SU(A') symmetry. To access the corresponding phase diagram experimentally, a set-up relying on atoms 
in optical lattices has been proposed recently [33]. These models enjoy a global U(A') symmetry, but the 
U(l) part of it is gauged with non-dynamical gauge fields, so that the Gauss constraint is the statement 
that the charge density, and therefore charge, vanishes identically. This means that a chemical potential 
for U(l) has no effect on the theory, as the U(l) charge is forbidden. We will see a manifestation of this 
in the dual variables as well. Generic chemical potentials can be rotated to those coupling to the — 1 
diagonal generators of the symmetry. The corresponding lattice action reads (by ^j,j = 1,2...N we 
denote the chemical potentials for the N components). 


S[z,U] 


N d 

—J X] X/ X/ ^"’‘^Zj{x)* Uy{x)zj{x + v) + ^''•‘^Zj{x) Uy{x)*Zj{x + v)* 

j=l xGA u=l 


(4) 
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The degrees of freedom are A^-component complex vectors with components Zj{x) € C. The spins are 
normalized to 1, i.e., = 1- As for the 0{N) model we consider a d-dimensiona \Nf-^xNt 

lattice A with periodic boundary conditions. The nearest neighbor terms of the action are connected 
with U(l) link variables, Uu{x) = with Ai,{x) G [0,27r]. In the temporal direction (z/ = d) the 

terms are non-trivial and give a different weight to forward and backward propagation in time. 

The partition sum is given by Z = f T>[C/]P[T] exp(—5[T, f/]) and the coupling J is positive (ferro¬ 
magnetic case). The integration measure is the product over the measures for the z{x) at each 
lattice site and 'D\JJ] is the product of Haar measures for the Ui,{x) G U(l) at all links (see (|2T]i and 
p2p below for the explicit form). 


3 Dual representation of 0{N) models with a chemical potential 

Both actions ([2]) and ([4]) obviously have a complex action at non-zero chemical potentials. In this and 
the next section we discuss the exact mapping of the corresponding partition sums to dual variables. We 
begin with the O(A^) model. Some of the steps here will then also be used in the next section for the 
CP(A^ —1) model. 

As the first step towards a dualization we introduce an explicit parametrization of the spin components 
ri 2 (x)e*‘^^^^ and f±{x) from ([2]) in terms of generalized spherical coordinates, 

^12(3^) = sin 0 Ar(a:)... sin04(a:) sin03(a:) , 

r3(x) = sin0Ar(x)... sin 04 (x) cos 03 (x) , 

r4(x) = sin0Ar(x)... sin 04 (x) cos 04 (x) , 


rN-i{x) = sm6j\f{x)cos6N-i{x) , 
XNix) = cos 9 j\i{x) . 


(5) 


The ranges of values for the angles are 4>{x) G [0,27r) and Oi{x) G [0, tt]. The path integral measure is a 
product measure over each individual spin and in terms of the generalized spherical coordinates has the 
explicit form 


Jv[f] = Y[J^^d^r{x), ( 6 ) 

sin03(x) (i03(a:) [ sin 9 ^{x)'^d64{x)... [ sin 9 ^{x)^d 9 n( x) . 

Jo Jo 


d^r(x) = 


IgN-l 

d(f>{x) 
1 2-k 


We remark that an irrelevant normalization constant was dropped. 

For a convenient notation the action is written as the sum of terms Sv\x) given by 

5i^)(x) = ^ ri2(x) ri2(x + u) , 

(a;) = ^ ,^^2 (a;) ri 2 (x + z)) g*(0W-0(x+P)) ^ 

S^^\x) = J rj{x)rj{x + v) , j = ‘i,A...N . 


(7) 


The action thus reads S[f] = —J2x,u,j Su^\x), where each term Sy’{x) is a simple product of neigh¬ 
boring spin components. 
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We now factorize the Boltzmann weight and expand each factor individually: 


Z = 


J Vlf] e = Jvlr] = J T^[r] 


= W]I[ 




E n 
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^ oo 
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^X,U 


u=0 
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iv .xj (si;\ 
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(f) I 


N J 

j "'X.l/ 


k ^ 1 » 
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xxO ) I 
Kj=3 ^x,i'‘ 




( 8 ) 


The individual exponentials of the terms Su"’ (x) of the action have been expanded with summation 
variables attached to the links of the lattice. We use kx^i^ € Mq and € INq for the first two terms 

[j = 1,2) and G Mq for the terms with j = 3,4... N. In the last line of ([8]) the explicit expressions 
for the (x) have been inserted. After ordering the various factors, the remaining integrals over the 
angles at the individual lattice sites are collected in I[{k,l,n^^^}]: 


\ Jo 27r Jq 


l+a. 


( 2 ) 


(9) 


(4) 

X d 64 {x) COS 94{x)°'^ sin04(x) 

Jo 


2+ai^4a® 


(N) ,r n , (2) , (3) (iV-1) 

d9Nix) COS 9 n{x)°'^ sin0Ar(x)^ , 


where we use the following shorthand notation 

= E + ^x-u,u + lx,u + lx-u,iy] , a-x^ = E + nx-u,u] ! J = 3,4 ... iV . (10) 


The integrals over the azimuthal angles (/)(x) give rise to Kronecker deltas which in turn enforce the 
constraints '^y[kx,v — lx,u — ikx-i>,u — lx-u,u)] = 0 Vx. The integrals over the polar angles are of the 
form (a, b G Mq) 


^ r(2^)r(l^) 

Iia,b)= rd0(cos0)“(sin0)^ = J r(2±^) 

Jo 1 0 


a even 
a odd 


( 11 ) 


(i) 

and therefore also impose constraints which restrict the admissible configurations of the variables rix/u G 
INq. The structure of the constraints can be simplified by the following change of variables: 

^x,u ^x,i/ — xfix^u , '^x,i' ^ ^ and kx^u T ^x,v — T 2 tHx^u 7 '^x,i' ^ 1^0 * (12) 

The partition function in dual representation assumes the form 


Z = 


E n^(E 


'^X,L' ^X — U,^ 


N 

n E (»“) 


(13) 


i=3 


(J/2)l"ia:.i^l+2 J rSJ], 


„ 7,^ i H 


{\mx,u\ +mx,y)\ rux.j Mi 


j=3 nx,u- J \ X j=3 


N j-l 

nn^(«-^.i'-2+E' 




k=2 
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with the a 


(j) 


now given by 


a 


( 2 ) 


^ ^ “h \'^x—v,v\ “h '^x—u,v)\ 

V 


,(i) 


I] [4^i + , j = 3,4 ... TV . 

(14) 


In its dual form the partition sum Z is a sum over the dual variables mx^v € Z, rux^v G INo 
f?) 

nx'v G INojj = 3,4...A^. These are subject to constraints at each site x. The constraints are written 
explicitly in the first line of (fT^ using Kronecker deltas b{n) and the evenness function E{n): 


6{n) 


1 n = 0 
0 else 


E{n) 


1 n even 

0 n odd 


(15) 


The second line in (fT311 gives the weight of a configuration {m,rn,n^^^: The first factor contains the 
coupling to the chemical potential, which now obviously is real and positive in the dual representation. 
The second factor in the weight comes from the expansion of the exponential functions and the third 
factor is the contribution from the integral over the local spin variables. 

Note that by explicitly writing the constraints in the first line of ([T^ the ^ are forced to be even 
by the evenness functions E[ax^) and we can omit the distinction of the two cases in (llip . using its 
first line only. Inserting these integrals in the product in the last term of ([T^ one finds that one term 
in the numerator cancels the denominator of the previous factor, leading to the final form of the dual 
partition function 


^ = E (n«(E 


N 

^x,i/ ^x—iy,u\^ E 

J=3 


(i)' 


(16) 




'^x.d 


-pj- (j/2)I"i-.-I+ 2 "i-.- 
y {\mx^u\+mx,y)\ mx,y\ n^J^] 



nf=. 

>ri 

(L 

^1 + bj^2 T 

)) 

ri 




where the ax^ are given in dMl). 

The final result (fl6]l provides an exact reformulation of the partition function Z in terms of the dual 
variables nix^u £ 'itt>x,u £ INo 3'3d nx,u € INoj j = 3,4... A^. All weight factors are real and positive for 
arbitrary values of the chemical potential fi and Monte Carlo simulations are possible without complex 
action problem. Large values of the dual variables are suppressed by the factorials coming from the 
expansion of the exponential function (and the T-factors do not change this). 

It is important to note that the configurations of the dual variables are subject to constraints for the 

(i) 

rux^u and the nx,u,j = 3,4... A^. For the variables nix^u the Kronecker deltas imply 


- m-s-i),!/] = 0 Vx ^rhx = 0 Vx . (17) 

V 

Obviously the Ihs. of that equation is a discrete version of the divergence Vm^, and the constraint simply 
implies that the discrete current rfix € is conserved. Thus the admissible configurations of the 
variables mx^v (he., the configurations that obey the constraints) are closed loops of flux. 

This also gives rise to an interesting geometrical interpretation of the terms coupling to the chemical 
potential: In (flbll the chemical potential couples to Yl,x''^x,d, i e-> the total flux of mx,u in direction 
V = d (i.e., the time direction with extension Nt). Since the m-flux is conserved the expression J2x''^x,d 
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corresponds to Nf times the total winding number of the loops of m-flux around the compactified time. 
Thus we find that the chemical potential /x couples to Yhx'^x.d = NtQ, where the net particle number 
Q is identified with the total winding number of the loops of m-flux around compactified time. In the 
bosonic case at hand this number is not restricted, but can take any integer value. Moreover, writing 
as the curl of another integer-valued field, one can show that the chemical potential couples to the 


m 




topological charge counting kinks in the dual Sine-Gordon model 

For the other set of constrained variables rix'u G = 3,4... A^, the constraints imply that the 

have to be even for all x and j = 3,4... A^. This can be written as 


sums Ox ^ 


x—u^ux 


H [4-',I+ <-£>,J even Vx 


,(i) 


U 


even Vx 


even Vx . (18) 


In the second step we have changed the plus in the individual terms to a minus which corresponds to 
adding an even number on both sides. As a result the constraints can again be written as a divergence 
and the constraint implies that the divergence Vax \ j = 3,4... N has to vanish only modulo 2 (which 
corresponds to a Z 2 symmetry). 

4 Dual representation of CP(A^—1) models with chemical potential 

Again the first step in the dualization is to write the spins z{x) G with \z{x)\'^ = E^=l \zj{x)\^ = 1 
using suitable coordinates. In particular we write the components Zj{x) with polar coordinates in the 
complex plane, 


Zj{x) = rj(x) with ^Pj{x) G [0,27r) , rj{x) G [0,1] and 


N 


Y^rjixf = 1 - 

i=i 


(19) 


The moduli rj{x) are subject to the normalization Xj{xY = 1 which restricts them to the N -1- 
sphere as in the 0{N) case. However, in addition we have the condition rj{x) > 0 in order to obtain a 
unique representation with the polar coordinates (fl^ chosen in the complex plane. We can parameterize 
the rj{x) again with generalized spherical coordinates, 

ri(a:) = sin0Ar(a:)... sin04(a:) sin03(x) sin02(3^) j 
r 2 {x) = sin 6 < 7 v(x)... sin 6 * 4 (x) sin 03 (x) cos 6 * 2 (x) , 
rsix) = sin 6 < 7 v(x)... sin 6 * 4 (x) cos 03 (x) , 


r]\f-i{x) = sin0]\f{x)cos9N-iix) 

rNix) = cos9n{x) , 


( 20 ) 


but now all angles 9j{x),j = 2,4... A^ run only over the interval [0, 7 r/ 2 ], such that the rj{x) are non¬ 
negative. With this parameterization we can write the path integral measure T)[z] explicitly as (we again 
omit parts of the normalization which drop out in expectation values anyway). 


'm =nl n 

V-i 



d9k{x) cos9kix) sm{9k{x)) 


2k-3 


( 21 ) 
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For the CP(A^—1) case we decompose the action (HD in the form 


S[z,U] 


N d 

j=l X u=l 

j ^’',dzj{x)* Uii{x)zj{x + v) , slf\x) 


( 22 ) 


Je>^i^-'’'^Zj{x) U^{x)*Zj{x + i>y . 


The link variables Uy{x) are for now considered as fixed external fields (we integrate them out later). 
For finding the partition sum Z[U] in a background configuration of the link variables we again factorize 
the Boltzmann weight and expand each term individually: 


Z[U] = lv[z] = lv[z] n = /P[T] E 




(?) 

(sW(x))'-'' » (s«w) 








uU) I 


E- 


k^i=o 


^X,V' 


U) 


E n 


J 


{fcO),P'T 


uU) I tO’) I 

i^XM' ri^ 




\X,U 


By -y)j we denote the sum over all configurations of the expansion indices ki^l,k^xl G INo 

assigned to the links of the lattice. The remaining integral over the powers of spin components at 
each site of the lattice can be factorized into two parts, k^^'^}] = k^^^}], 

which correspond to the integration over the phases and the moduli of the spin components. 

The integral over the phases ^Pj{x) is rather trivial: 




n 


1 27r 




U) 


-uU) 'll 

x — i>,vk\ 



(24) 


it simply gives rise to a product of Kronecker deltas for all sites. At every site the Kronecker deltas 
enforce the constraint 


=0 , i = l,2...7V. 


(25) 


The integral over the parameterized moduli rj{x) is given by 


z,[{fc«t“>}]=nn 


X k=2 


'0 


d9k{x) COS 6k{xY^°'^^ ^ sin6k{x)^^ ^+Sj=i 


(26) 


where we use the abbreviation 


^ ~ E \kxl> + + k^x-v,u + k, 


rSi) 

X — VM 


j = l,2...iV. 


(27) 


All the integrals that appear in f|26D are related to (ITTIl but the integration runs only up to 7r/2. Thus 
the result is always given by the first line of (llip with an extra factor of 1/2. When taking the product 










of all the integrals that build up ()26l) . again the numerator of a term and one of the denominators of the 
subsequent term cancel. The integral assumes the simple form (we dropped an irrelevant 

overall factor 2“^^) 


= n-^ 


nf.ir(^ (2 + oS?>)) 


(28) 


r(i(2JV + E"i<.S?’)) ' 

Note that here the integrals over the 9j{x) did not generate any additional constraints - the constraints 
(1251) from integrating the ^Pj{x) already imply the evenness of the o!i\ 

As before we can simplify the expressions further by rewriting the dual variables k!il € Mo in 
terms of new dual variables G Z and mx}y G Mq, 


^ and k^^), + k^^l = + 2 

The sums defined in (|27]1 turn into 


^ [\^xl\ + + k: 


I + 


j = 1,2...N , 


(29) 


(30) 


and we find for the partition sum Z[U] in the background field of the Uy{x)'. 

+ mxl)’ nlxll 


m = E TldE 




n 






n'i.r ( 1 ( 2 +„<?•))) 


(1 (2Ar + E''.,<iS?*)) 




E, 


‘'X.U 


\X,U 




,X,J 


(31) 


The final step0 is to integrate the link variables Uy{x) = with the measure J'D[U] = Hx ^ Jq^ 

which gives rise to another set of constraints: 


-T* 7/ -r TV *^0 ^ 


Aj) 


= nME 


m. 


U) 

X,V 


(32) 


The final result for the CP(A^—1) partition sum in the dual representation thus reads 


(i) _ ] 

X,U X — 0Mi 




nME 


m. 


U) 

x,v 






jj_nf=ir(^ (2 + a«)) 


with the aA given by (|3Ql) . Only the dual variables mx,l G Z are subject to constraints, while the 
m'J}, G Mq are unconstrained. All constraints are collected in the first line of (|33ll : The first factor 


(j) 


(33) 


,j i\mxl\ + mx,l)l mx}v\ j \x T (2N + a* ^)) / 


^Actually this could have been done at any stage of the derivation. 
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implies that all discrete currents are divergence-free, i.e., Vfhx'^ = ~ ''^x-u u\ ~ ^ 

similar to the 0(N) case in Eq. (I17|) . We stress that the condition of vanishing divergence holds 
individually for all j = 1,2... N. A second constraint ties together the components j = l,2...N: At 
each link the sum '^jrni^u of the dual variables mi^u G ^ has to vanish. Note that this implies that 
for the case of identical chemical potentials = //2 = = A model is independent of /i. 

This is just reflecting the fact that the chemical potential for the U(l) gauge symmetry cannot affect 
the theory. The reason for this is that the gauge field is non-dynamical, and the Gauss law is simply 
a statement that the charge density of this U(l) symmetry vanishes at every point, and the total U(l) 
charge must be zero for all physical states of the model. 

The weight factors are collected in the second line of psp . It is obvious that all of them are real and 
positive also for arbitrary values of the chemical potential and the complex action problem is solved in 
the dual form. As for the 0{N) case, the chemical potentials couple to the time component of the 
currents and combined with the flux conservation (i.e., the condition Vrux'^ = 0) this implies that the 
charge Qj the chemical potential fij couples to is given by the winding number of the m^-^^-flux around 
the compact time direction. 

We remark that in two dimensions {d = 2) one may consider a topological term by including an 
additional Boltzmann factor with Q[U] = {2Tr)~^J2x^'‘^Ui{x)U2{x + i)Ui{x + 2)*U2{x)*. This 

introduces an additional source for an imaginary part in the action. It turns out, that this latter complex 
action problem cannot be removed completely in our dualization (the same conclusion was found also 
in the alternative dualization in [H]). This is an interesting finding, since for U(l) lattice field theories 
that contain the same topological term but also an action term for the gauge fields, the dualization does 
indeed solve the complex action problem from the vacuum term [341135] . 

5 Remark on the relation of the dual 0(3) and CP(1) models 

An interesting question is how the dual representations of the 0(3) and the CP(1) models are related. 
These are the lowest-A^ cases of the two models and in the conventional continuum representation can 
be transformed into each other (in this discussion we restrict ourselves to vanishing chemical potential). 
In the continuum one can show via r/ = z*{ai)ijZj with a the Pauli matrices, that the 0(3) action 
(9i/f)^/2 turns into twice the quartic CP(1) action {d^z)'^{d^z) + (z'^d^z)'^. The latter occurs from the 
action with gauge fields {D^z)^D^z), Dy = + iAy, after integrating out the gauge field. 

On the lattice one can do the same. Under the field mapping above, the lattice 0(3) actions ([I]), 
(|2]1 turn into actions quartic in the CP(1) fields closely resembling the continuum quartic CP(1) action. 
On the other hand, the Boltzmann factor with the CP(1) action ([4]) (set N = 2), after integrating out 
the link field can be shown to be a product of Bessel functions and not of the form exp(—5) with S 
a polynomial action in the CP(1) fields. Hence, the 0(3) actions from which the dualizations start are 
different and the same is expected for the dualized partition functions. The differences should of course 
vanish in the continuum limit. Note that both dualizations possess three dual variables, (m, m, in 
the 0(3) case vs. in the CP(1) case (having solved the second constraint by putting 

These two dualizations for the 0(3) model could have different computational features 
and might highlight different aspects of our understanding of dualizations in general. 
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6 Summary and comments 


In this paper we have derived dual representations for the 0{N) and CP(A^ —1) models with chemical 
potentials. The dual representations use N sets (respectively 2N sets for the CP(A^ —1) case) of 
integer valued dual link variables. Some of the dual link variables obey constraints, in particular they 
have vanishing divergence, i.e., conserved flux. The chemical potential couples to charges which in the 
dual picture can be interpreted as the winding numbers of the corresponding conserved flux around 
the compact time direction. All weight factors are real and positive and in the dual representation the 
complex action problem is solved. Suitable dual Monte Carlo strategies were discussed in the literature 
(see, e.g., (THU, pHE] for examples) and for the case discussed here we are currently preparing a 
paper [35] presenting the results of a dual simulation in the representation derived here. 

Vacuum expectation values of observables can be obtained easily in the dual representation: Deriva¬ 
tives of InZ with respect to the couplings J and fj,j give rise to bulk observables and their moments 
which can be used for determining the various transitions in the models [36]. In addition one can go 
through the same dualization steps also for chemical potentials that depend on the space-time 

arguments x. These then may be used as source terms and after taking suitable derivatives of InZ one 
obtains the corresponding n-point functions. 

In the introduction we have discussed that there is no general strategy for finding a dual repre¬ 
sentation. Studying various dual forms of different models thus provides new tools for attacking more 
challenging non-abelian symmetries, in particular gauge theories with non-abelian gauge groups. How¬ 
ever, not only finding more tools for dualization is important, but also understanding the dual variables 
is an interesting issue per se. Here we have demonstrated that different dual representations exist for the 
same model, (in detail we discussed the case 0(3) = CP(1)), and one thus might try to analyze what 
constitutes the key ingredients of a dual representation. Understanding these key characteristics might 
open the possibility of formulating lattice field theories directly in the dual framework of matter fluxes 
and surfaces for gauge fields. Clearly this is a rather speculative perspective at the moment, but it might 
become more relevant with every successful dualization of a lattice field theory. 
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work is partly supported by the Austrian Science Fund FWF Grant. Nr. I 1452-N27 and by the DFG 
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